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Abstract: We show how to define a consistent braneworld cosmology in a model in which 
the brane is constructed as a field-theoretic domain wall of finite thickness. The Friedmann, 
Robertson- Walker metric is recovered in the region of the brane, but, remarkably, with scale 
factor that depends on particle energy and on particle species, constituting a breakdown of 
the weak equivalence principle on sufficiently small scales. This unusual effect comes from 
the extended nature of particles confined to a domain- wall brane, and the fact that they feel 
an "average" of the bulk spacetime. We demonstrate how to recover the standard results 
of brane cosmology in the infinitely-thin brane limit, and comment on how our results 
have the potential to place bounds on parameters such as the thickness of domain-wall 
braneworlds. 



Keywords: Large Extra Dimensions, p-branes, Cosmology of Theories beyond the SM 



Contents 

[l]. Introduction |] 

|2|. Fundamental brane cosmology ||| 

|3|. The extension to a domain-wall brane ^ 

|3.1| A localised scalar field 

3.2 Localised fermions |TT| 



|3.3| The effective Newton's constant 13 

HL Effective scale factor for a thin domain wall 15 

^ Conclusions 17 

|A] . Einstein tensor components 19 



1. Introduction 

Extra-dimensional theories of particle physics beyond the standard model have, in recent 
years, become a standard part of the array of phenomenological models that we hope to 
test with the imminent operation of the Large Hadron Collider (LHC). Some models are 
extended versions of the original Kaluza-Klein (KK) construction, in which the known 
matter fields propagate in a higher dimensional space, and our usual low-energy physics 
is augmented by the properties and interactions of their KK towers. More dramatically, 
inspired by the D-branes of string theory, Arkani-Hamed, Dimopoulos and Dvali (see 
also [§, and Randall and Sundrum (RS) ]|, |j| have suggested an alternative set of 
extra dimensional constructions, in which a restricted subclass of fields (often just gravity) 
propagate in the additional directions, while standard model particles remain confined to 
a 3 + 1 submanifold, or brane (see also || f7|). 

In such models, the brane is typically modelled as infinitely thin - a delta-distribution 
source - and one then constructs a four-dimensional effective field theory of the standard 
model particles plus those propagating in the bulk, with their associated KK partners. A 
rich phenomenology results, with details depending on the particular type of extra dimen- 
sional construction employed (see || |9|, [H], [llj for alternate constructions, and [12, 13] 



for general overviews.) Predictions include a spectrum of KK particles expected to be 
produced at colliders and modifications of the Newtonian inverse square law at sufficiently 
small distances. Details regarding cosmological evolution were worked out in the initial 

series of papers m, n, m, m. m, m, where & 

was shown that there are deviations from 



the standard cosmology at temperatures above a particular normalcy temperature. 
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Over the decade following these discoveries, an immense body of work has been per- 
formed to tease out the particle physics and cosmological signatures of these models. Cos- 
mology in particular has received much attention [2C], with a great deal of the work focused 
on inflation, gravity perturbations and similar cosmological consequences of braneworlds 
which would be evident today. Attention has also been given to more exotic phenomena, 
such as colliding branes 21, 22, 24, 25 1, which is a unique aspect of brane cosmology that 



can have drastic effects on the very early universe. Upcoming precision experiments and 
missions in these fields, along with particle physics results from the LHC, promise increas- 
ingly stringent tests in the near future. Such precision allows us to probe ever more fine 
detail of extra dimensional theories, requiring an examination of the sensitivity of the tests 
to the specifics of the underlying construction. One such detail is the delta-distribution 
nature of the brane itself. 

It is certainly possible that the 3-brane on which we find ourselves is an actual D-brane 
(or orientifold plane) , and that our universe is embedded in a string-theoretic realisation of, 
for example, the Randall- Sundrum braneworld. In that case, it would seem sensible to treat 
the brane as infinitely thin for all purposes of low-energy particle physics and cosmology. If 
it is necessary to go beyond this approximation and model the effects of a finite-thickness 
fundamental brane, then one can do so by averaging the relevant cosmological quantities 
over the extent of the brane, and imposing boundary matching conditions where the thick 
brane meets the bulk ||, ||, |f, H- 

However, there is another possibility, that the braneworld model itself is unrelated 
to string theory, but is instead a field-theoretic construction, in which the brane is a 
topological defect (a domain wall, for example) in a higher dimensional space [30, EHl. 
Generally, a scalar field is responsible for generating the topological defect, and the brane 
and Randall-Sundrum warped metric become smoothed out versions of their counterparts 



in the fundamental case f32| , 33, [34|, 35, p6[ . In such scenarios, the extended nature of the 
brane in the extra dimensions may have important implications for phenomenology and 
cosmology, both through its complicated relationship to the bulk, and through the field 
theory mechanisms through which it must confine the standard model particles. 

In a series of papers, two of us (DPG and RRV) have examined this question in detail 
for the case in which the brane is a codimension one soliton (a domain wall) in 4 + 1 
spacetime dimensions, providing a Minkowski background for standard model fields. In 



such a construction |37| (see also [38] for an extension to a larger grand-unified group), a 
subtle interplay between the bulk gauge fields (and gauge group), the bulk scalar fields, and 
the symmetry breaking field that forms the domain wall is necessary for standard model 
fermions, the Higgs field, and the SU(3) x SU(2) x U(l) gauge bosons to confine and appear 
in the resulting 3 + 1 dimensional effective theory. 

In this paper, our goal is to extend previous work to understand the cosmology of 
field-theoretic braneworld models, requiring us to abandon the assumption of Minkowski 
spacetime on the brane, and to develop a framework in which the standard model fields 
feel a background Friedmann, Robertson- Walker (FRW) metric. Related work on thick 
fundamental brane cosmology [27, 29 1, and the prescription for averaging 5d quantities 
to produce the corresponding effective 4d quantities, is drastically modified in the case of a 
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domain-wall brane. As we shall see, the KK expansion of 5d fields and the integration over 
the full extent of the extra dimension leads to some interesting and unexpected features. 

The structure of this paper is as follows. In the next section we review the main points 
of brane cosmology as constructed around a delta-distribution source. In Sec. || we then 
demonstrate how to build the brane from a field-theoretic domain wall, and discuss how 
to define the metric on such an object. We then show how scalar fields and fermions are 
confined to such a domain wall, and derive expressions for the cosmological metrics they 
experience. The thin wall limit is verified in Sec. || and we then collect and discuss our 
main results. 



2. Fundamental brane cosmology 

In this section we summarise previous results for the cosmological evolution of a funda- 



mental brane with localised sources (see [14, 15, 18] for details, and ||39|| for an extension 



to higher codimension.) The idea is to take a 5-dimensional bulk spacetime, include brane, 
bulk and brane-localised stress-energy sources, and solve the 5-dimensional Einstein equa- 
tions. The brane is considered to be a fundamental object and is modelled by a delta 
distribution, with the total action being 

S = j d 5 X^ M$(R- 2A) + J d 5 2V-Sbrane <%)£brane , (2.1) 

where M5 is the 5-dimensional gravitational mass scale, A is the bulk cosmological constant 
and g and <?brane are the- determinants of the metric in the bulk and on the brane respec- 
tively. The delta-function localises £brane to the brane, which includes the brane tension 
and standard model fields. 

Since we are interested in cosmological solutions, we consider sources that are homoge- 
neous and isotropic in the three spatial dimensions, and the most general metric consistent 
with these symmetries 

ds\ = -n 2 (t, y)dt 2 + a 2 (t, y)j ij dx i dx i + b 2 (t, y)dy 2 , (2.2) 

where i,j run over the three spatial dimensions, and 7^ is the metric of the three-space, 
which may be positively curved, flat or negatively curved. The Einstein equation is 

1 1 

Rmn — -^9mnR + <7matA = 2 jy^3 Tmn i (2-3) 

where M, N are 5-dimensional indices. For explicit expressions of the components of the 
Einstein tensor, see Appendix [A]. The brane tension and brane-localised sources are rep- 
resented jointly by the 4-dimensional density p b and pressure p b and appear in the stress- 
energy tensor as 

T M N = 6(by)dmg(-p b ,p b ,p b ,p b ,Q) . (2.4) 

For general values of A, p b and p b , Einstein's equations will yield time-dependent 
solutions, corresponding, for example, to an expanding spacetime on the brane. Before 
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exploring such solutions, we note that it is possible to fine tune the sources to produce 
a time-independent background. This corresponds exactly to the scenario of Randall and 
Sundrum [||, |B| who demonstrated that gravity is localised to a fundamental brane tuned 
in such a way. 

The specific choice necessary is that the brane source be pure tension (corresponding 
to a 4d cosmological constant) pb = —pb = cr, tuned against the bulk cosmological constant 
according to 



-24MfA . (2.5) 

Note that this implies that the bulk geometry must be five-dimensional Anti-de Sitter space 
AdSs, with A < 0. The corresponding metric solution is then 

dsl = e- 2M (-dt 2 + jijdxUx^ + dy 2 , (2.6) 

where 

/i= V? = l2^f (2 ' 7) 
Moving back to the general time-dependent case, one solves the 5-dimensional Einstein 
equations by imposing the Israel matching conditions - calculating the discontinuities in 
the derivatives of the metric components across y = 0, and relating these to the delta- 
distribution sources (see |l4T| ). It turns out that the behaviours of pb, Pb and the metric 
components evaluated at y = are independent of the metric solutions in the bulk, and 
obey 

p b + 3H { Pb + p b ) = 0, (2.8) 

F ° 2 = l4&| + i-| + 4' (2 ' 9) 

where the parameter k takes values +1, 0, —1 according to whether the metric 7^ describes 
positively-curved, flat, or negatively-curved spatial 3-sections. Here a dot corresponds to 
a time derivative and time has been rescaled such that no(t) = n(t,y = 0) = 1. The 
integration constant C represents an effective radiation term, or so called "dark radiation" 
(see |l8| for bounds on this term from nucleosynthesis) and from now on we will set C = 0. 

The effective 4-dimensional scale factor ao is the 5-dimensional metric component 
a(t,y) evaluated on the brane ao(t) = a(t,y = 0) and Hq is the corresponding Hubble 
parameter. Equation ( [2.8| ) describes the usual 4-dimensional conservation of energy on the 
brane (the continuity equation) and Eq. ( |2.9|) is the modified Friedmann equation. Due 
to the proportionality Hq ~ pb instead of the usual Hq ~ y/pi, this Friedmann equation 
seems at odds with observation. The clue to fixing this problem comes from considering 
the time-independent Randall- Sundrum solution, where the brane tension contributed an 
energy that exactly cancelled the bulk cosmological constant. Guided by this, one writes 
the total brane source pb, Pb as a sum of a background brane tension a and some other 
general brane source p, p 

p b = a + p, (2.10) 
p b = -a+p, (2.11) 
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where a is defined as in fl2,5|) . The effective Friedmann equation for ao now reads 

Hi = -Z-zp + —^p 2 - 4 • (2-12) 
72Mf F 144Mf P a 2 v ; 

If we assume p is small compared to a, then the p 2 term gives small corrections to the usual 
behaviour and the evolution of ao is driven to first order by p, with the proportionality 
constant playing the role of the effective Planck mass via 

, o 12M| 6M. 3 

Ml = b - = 5 . (2.13 

An important feature of cosmology in codimension-1 braneworlds is that this entire 
analysis is independent of the behaviour of the metric components in the bulk. Nevertheless, 
it is possible to find bulk solutions, and since we will make use of them in a later section 
we provide them here. For C = they read (l^, [13| 

n(t,y) = - esinh(p\y\) , 

a(t,y) = ao(t)[e-^\ - esmh(p\y\)] , (2.14) 
b(t,y) = 1 , 



with p defined as in (|2.7|) and 



e^P, (2.15) 
a 

e = e + -^. (2.16) 

Note that for e = and k = we recover the RS warped-metric solution given by Eq. (^^) . 

The parameter e measures the energy density in matter and radiation, relative to the 
tension of the brane. In terms of this parameter, the Friedmann equation ( |2.12| ) is 

demonstrating that e 1 is required to recover conventional cosmology. The earliest 
direct evidence for our standard cosmological evolution is provided by primordial (big bang) 
nucleosynthesis (BBN), which takes place at temperatures of order an MeV. Therefore, we 
are safe from cosmological constraints if we choose a 3> (lMeV) 4 . 

Finally, we briefly discuss the extension of these results to fundamental branes with 



finite thickness [27, |2^, pS|] . I* 1 these scenarios the brane and brane-localised sources are 
modelled as stress-energies distributed over the finite thickness of the brane. The effective 
4d quantities, such as the scale factor, energy density and pressure, are defined to be the 
spatial average, over the extent of the brane in the extra dimension, of their corresponding 
5d quantities. One then rewrites Einstein's equations in terms of these averaged quantities 
and identifies corrections to the infinitely-thin brane scenario. This averaging prescription 
is an important first step in understanding cosmology away from the infinitely-thin brane 
limit. However, a more complete treatment is essential, since, for example, in the Minkowski 
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domain- wall set-up, one needs to expand the 5d fields in KK modes and integrate over the 
full extent of the extra dimension. The rest of this paper is devoted to the development of 
a more complete averaging framework, within which it is possible to analyse the cosmology 
of domain-wall brane scenarios. 



3. The extension to a domain-wall brane 

Our main goal in this paper is to extend the analysis of the previous section to the case 
in which the brane is topological defect - a domain wall generated by a scalar field. The 
central problem is how to identify the effective 4-dimensional scale factor (the analogue of 
ao) and the equations that describe its time evolution. As we shall see, this question turns 
out to have an interesting and nontrivial resolution, which may have specific implications 
for the signatures of such field-theoretic braneworlds. 

The creation of a domain-wall brane coupled to gravity is quite straightforward; we 
will follow closely the construction in |^(|. Beginning with a scalar field x an d a suitable 
potential, boundary conditions are chosen so that x develops a kink-like profile, which can 
be thought of as a y-dependent vacuum expectation value. As y — > ±00, the value of x 
approaches vacuum and its energy density rapidly approaches zero. However, due to the 
topology of the vacuum (in general a discrete symmetry is required), a domain wall forms 
around y = 0. The combination of gradient and vacuum energy in the core of this object 
plays an analogous role to the brane tension a in the fundamental case of the previous 
section. The shape of the distribution of stress-energy due to the y-dependent profile of x 
is a smooth version of the fundamental delta-function brane. In the non-cosmological case, 
i.e. when one seeks the Minkowski metric on the brane, the solution for the metric then 
yields a correspondingly smooth version of the e _M ' y ' warp factor in (|2.6| ). 

Because this domain- wall brane is extended in the extra dimension y, any fields that 
were previously brane- localised by the delta function are no longer strictly located at y = 0. 
Rather, such fields (typically the standard model fields) must first be written as full 5d 
fields, which are coupled to x 111 such a way as to produce a Kaluza-Klein tower of 4d 



fields on the domain wall (in the Minkowski-brane example see e.g. [41 1). The ground 
state profile of the 4d tower has a Gaussian like shape which, when squared 1 , reduces to a 
delta-distribution in the limit of an infinitely-thin domain wall. 

Here we are interested in the more general cosmological case. Our objective is to 
understand the effective 4-dimensional metric on such a domain-wall brane and how the 
localised fields propagate in that spacetime. In the fundamental-brane case, 4d fields are 
located at exactly y = and have no y degrees of freedom. The 4d metric they feel is 
thus just the 5d metric evaluated at y = 0, and for the RS metric ( |2.6[) this slice is just 4d 
Minkowski spacetime. For the cosmological metric fl2.2|) the slice at y = has the form 

ds 2 = -n 2 (t,y = 0)dt 2 + a 2 {t,y = 0)j ij dx i dx j . (3.1) 



1 The squaring of the extra-dimensional profile comes from the normalisation of the kinetic term, which 
is quadratic in the field, hence quadratic in the profile. 
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By scaling t such that n(t,y = 0) = 1, it is clear that the effective 4d metric is of the 
FRW form, with the effective scale factor defined by a e f[(t) = ao(t) = a(t,y = 0). The 
solutions to the 5d Einstein equations given in the previous section then describe how a e g 
evolves, and hence describe the spacetime in which the localised fields propagate. In this 
fundamental-brane scenario, each field has the same time (with the same normalisation) 
and feels the same scale factor, and so it is sensible to say that the effective 4d metric is 
unique and defined by (|3.1|). 

For the domain- wall brane scenario things are quite different and, as we demonstrate 
explicitly below, we are led to abandon the question "what is the effective scale factor on 
the brane?", and allow that different fields may propagate in different spacetimes. The 
essential reason for this comes from the extended nature of the profiles, as the associated 
fields are now sensitive to the metric around y = 0, not just the slice exactly at y = 0. Since 
the cosmological evolution of the slices in the vicinity of the brane are misaligned (they 
expand at different rates), there is a kind of "dimensional parallax" effect, whereby different 
species of particle are subject to a different averaging (they have a different perspective) 
of the slices. 

We note here that for the Minkowski domain-wall brane with a smoothed-out version 
of the RS metric (|2.6| ), things are much simpler, because each 4d slice is proportional 
to Minkowski spacetime. Therefore, the Minkowski part essentially factorises out of the 
averaging integral and each field feels the same spacetime. 

The analysis in Sec. || determined the effective scale factor ao in the case of a general 
brane-localised source parameterised by p and p. For the domain-wall brane scenario, we 
need to look at the sources from the more fundamental level of classical fields. The general 
strategy is to identify the kinetic term in the action for the relevant field, integrate out the 
extra dimension y, and then to match the resulting 4d effective action to the canonical 4d 
action for such a field in an FRW background. 

3.1 A localised scalar field 

We begin by considering a scalar field, turning to fermions in the next subsection. We take 
the 5d metric given by Eq. ( [2. 2D and a 5d scalar field <fr(t,x l ,y) separated, for reasons we 
shall expand on below, as <I>(t,x*,y) = f(t,y)(j)(t,x l ). The objective is to determine the 
effective 4d spacetime on which the relevant 4d field <f> propagates. 
The action for a 5d scalar Q(t,x l ,y) with metric gMN is 



~g MN d M $d N <S>-U(<S>) 



(3.2) 



where the potential U may contain couplings of <3? to the domain-wall field (to localise <&) 
or couplings to other fields. Using the metric ansatz (|2.2| ) we then obtain 

' . (3.3) 



J d A x J dy na 3 b^ ~( - n' 2 ® 2 + a' 2 ^' d,$ + b~ 2 {d y <$>) 2 ) - U($) 



By analogy with the flat case, our first instinct might be to separate variables by writ- 
ing <&(t,x l ,y) = J2 n fn(t, y)-7r n (x l ). However, in the case of a time-dependent metric such 
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an expansion makes it difficult to identify a 4d scalar field, since the time dependence of, for 
example, a 4d plane wave, is consumed by the profile f n , and ir becomes merely a static spa- 
tial wave. The next obvious suggestion is to instead write $>(£, x l ,y) = Y^ n fn(y)4>n{t, x l ), 
so that (j) n can be identified as a 4d Kaluza-Klein mode with extra dimensional profile f n . 
Here however, we encounter a different problem, namely that the time variation of the 
metric components implies that the extra-dimensional profile will in general change with 
time. 

We overcome these obstacles by noting that there are really two time scales in the 
problem: the cosmological time scale of the evolution of the background metric, and the 
time scale associated with particle physics processes. With this in mind, we consider the 
following separation of variables 

1>(t,x i ,y) = Y,fn(t,y)Mt,x l ) . (3.4) 

n 

The possible ambiguity in the time dependence (whether it appears in f n or <f> n ) is resolved 
by the requirement that 4> n satisfies the 4d Euler-Lagrange equation, which will be specified 
below. In order for cp n to be identified as a propagating 4d field, it must also carry the 
majority of the time dependence, hence we impose the condition f n /f n <C 4>nl4>n- These 
requirements formally identify the class of solutions for f n that we are allowing. 

One should consider this prescription a separation of scales, rather than a strict sepa- 
ration of variables. Quantitatively, fin/fin ~ E where E is the energy of the particle, and 
fn/fn ~ H where H is the Hubble constant. In natural units we have H ~ 10 _32 eV, which 
is tiny compared to the typical energy of a particle. In what follows, we therefore neglect 
all time derivatives of f n and of the metric components n, a and b, since they are much 
smaller than the other terms in the action. 

From now on we focus on a single mode of the KK tower and drop the subscript n. 
Then, with the prescription ( |3.4[) and the assumptions regarding small time-derivatives, 
the action becomes 

S 5 = y>V7 |%[-^(-^/ 2 ^ + na6/VW^+^/V)-na 3 6C/] , (3.5) 

where a prime denotes a derivative with respect to y. The third term, proportional to 
(j) 2 , will contribute to the potential U. Integrating over the extra dimension yields the 4d 
effective action 



±(-F(t)fi + G(t)jVd i( f>d j( />) +... 



(3.6) 



where we have written only the kinetic terms explicitly, and defined 

F(t)^Jf 2 ^dy, G(t) = J fnabdy. (3.7) 

The action (3.6) is almost what we are looking for, but what remains is to correctly 
identify the 4d line element describing the spacetime within which (j) propagates. To do 
this, we match to the prototype line element 

ds\ = -T 2 (t)dt 2 + X 2 (t)- /ij dx i dx j , (3.8) 
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and the corresponding prototype action 

^(proto) = J d * xT{t)X 3 {t) ^ _ 1 (_r~ 2 (t)0 2 + X-^t^di^) . (3.9) 

Matching the effective action ([O]) with the 4d prototype ([O]) we then obtain 

F(t) = T _1 (t)X 3 (t) , G(t) = T(t)X(t) . (3.10) 
Solving for T(i) and -X"(t) gives 

-1/4 / „ s 3/4 

/ / _o \ 



r(t) = F-Ww&'Ht) = (| V)~ (/ AaWl/) , (3-11) 

X{t) = FV\t)GV\t) = (| / 2 ^*/) V4 (/ fnabdyX'' . (3.12) 



The time-dependent functions T(t) and X(t) define, along with (j3~lj|), the effective 4d 
line element followed by the field eft. As we shall soon demonstrate, we are free to rescale 
/ by an arbitrary (slowly varying) function of time, and we can use this freedom to fix 
T = 1. This corresponds to choosing a canonical time coordinate. The scale factor for eft is 
then precisely 

H (t) = X(t) . (3.13) 

Notice that the temporal behaviour of X(t) (and hence cu) is inherited from the time- 
dependence of the metric components and possibly f(t,y), all of which are taken to be 
slowly varying. 

This result for the effective scale factor immediately raises two important points. The 
first is that scalar modes with different profile functions will have different definitions of the 
scale factor a^. Thus, it is not possible to define a unique scale factor for this 4d effective 
theory. Instead, each 4d scalar field, whether it arises from a different 5d field, or is merely 
a different KK mode of the same 5d field, propagates according to a different effective 4d 
metric. 

The second interesting point is that the procedure above will yield a different result 
for a fermionic field (and also other spin fields) due to the difference arising from the spin 
connection in the kinetic term. We will follow this point up in the next section where we 
explicitly perform the relevant calculation for a fermion. 

As a consistency check, we consider Eqs. (3.11) and ( p. 12 ) in the limit of an infinitely- 



thin domain wall. In such a limit, the square of a typical ground state profile / becomes 
proportional to a delta- function distribution, f 2 — > S(by). This comes from the kinetic 
term for 0, which is quadratic in /, and must be normalised such that, in the thin brane 
limit, / f 2 d(by) = 1. The integrals for T(t) and X(t) can then be performed analytically 
yielding T{t) = n(t,y = 0) and X(t) = a(t,y = 0). These coincide with the fundamental 
brane case, in which the 4d metric is the 5d metric evaluated on the brane. 

A further check on our derivation can be made by looking at the separable (but less 
general) version of the cosmological metric, given by 

ds 2 b = c 2 (y)[-dt 2 + o^thij-d&W] + b 2 {y)dy 2 . (3.14) 
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This ansatz allows for A0IS4 and CIS4 brane solutions, as detailed in [|32|, |42|, 43, |44j, |45|, 
||, @. The effective 4d metric for <t> then has T(t) = (// 2 c 2 S (fy) 1 / 2 and = 
a(t)(J f 2 c 2 b dy) 1 / 2 . Note that T is constant (/ will be time- independent; see later) and 
we can normalise / to make T = 1, and then find that X(t) = a(t). In this case we again 
recover the known result, namely that all fields on the brane feel the same metric. 

To complete this formal analysis of we determine the differential equation satisfied 
by the profile function f(t,y). Our definitions above for the separation of scales ensure 
that behaves as a 4d scalar field in a spacetime characterised by T(t) and X(t). This 
means that (f> will satisfy the Euler-Lagrange equation 



+ W " T W :&<!>) = m 2 ^ , (3.15) 



where r .. are the connection coefficients associated with the 3-space metric 7™ and m is 
the effective 4d mass of 4>. The parenthesised term on the left hand side is simply the double 
covariant-derivative of 4> with respect to 7^. Note that we are ignoring time derivatives of 
T(t) and X(t), which are much smaller that the derivatives of (f>. 
Now consider the 5d Euler-Lagrange equation for <E> 

dU 



g MN (d M d N * - r (5) p MN d P <s> ) = ^ , (3.16) 



(5) p 

where T mn are ^ ne ^ connec ti° n coefficients. We first separate variables, neglect time 



derivatives of n, a, b and /, and use Eq. ( 3.15| ) to eliminate the spatial derivatives of 



(thus the 4d mass will appear). We then linearise the equation, yielding 



where dU/d<f> = + C($ 2 ). Notice the appearance of the <f> term, which is absent 

when we specialise to Minkowski spacetime on the brane. There are two reasons for this. 
First, there is a mismatch between the 5d ratio of the time and 3-space metric factors, 
and the corresponding effective 4d ratio; a 2 /n 2 7^ X 2 /T 2 . For a Minkowski brane these 
ratios are equal because each 4d slice of the 5d metric is proportional to Minkowski space- 
time. Second, we have employed a separation of scales rather than the usual separation 
of variables (which did not work in this setting). This <p term then quantifies the inability 
of the domain-wall brane to localise proper 4d effective fields, at least in a cosmological 
background. 

To proceed, we need to eliminate the <p and <j> factors so that we have an equation that 
can, at least in principle, be used to solve for /. To this end we solve the 4d Euler-Lagrange 
equation ( 3.15| ) (in flat space; k = 0) and find "plane waves" of the form 

(j)(t, x { ) oc exp(-iT 2 Et + /.V 2 - , ; //.r ; ; , (3.18) 

where E is the energy of the wave, p % is its momentum, and T 2 E 2 = X 2 ^Hjp L p > +m 2 . Then 
(j) = —E 2 T 4 (j), and equation (|3.17| ) becomes 

,X 2 „ m E 2 T 4 fX 2 a 2 - ~ 



f' + [ n - + "-- V r)f + b 2 

n a b 



m-^r - £/ (1) 



2_ 

a 2 a 2 \ T 2 n 2 



f = . (3.19) 
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Usually, such an equation depends only on m, implying that although different masses in 
the KK tower of 4d fields (cf>Q, (f>i, etc.) have different profiles, these profiles are independent 
of the energy. Here however, the equation also depends on E, so that quanta with the same 
mass but different energies (or momenta) have different profiles. On the surface, Eq. ( 3.1SQ 



looks linear and homogeneous in /, but it is in fact a non-linear integro-differential equation, 
since both T(t) and X(t) are defined in terms of /. Nevertheless, this equation still has 
the property that / can be rescaled by a y-independent factor, so long as the eigenvalues 
m and E are also appropriately rescaled to compensate for the change in T(t) and X(t). 
In fact, since / <C E, as discussed previously, we may even take this factor to have a 
(mild) time-dependence. As we advertised earlier, the rescaling of / can be used to choose 
a canonical time coordinate, corresponding to fixing T = 1, which is achievable precisely 
because T(t) depends on /. 

As a check on our derivation, for the separable cosmological metric fl3.14| ), the factor 
X 2 /T 2 — a 2 /n 2 vanishes, and Eq. (3.19) simplifies to the known result (see fl40|| ) 



m 2 \ - £/« 
c z 



f = . (3.20) 



Note the lack of time dependence, implying that / is a function of y only. The profile also 
no longer depends on the energy of the mode, just its mass, as usual. 

Let us summarise our results for scalar fields. Given a 5d background metric, described 
by functions n(t,y), a(t,y), b(t,y), and a generic coupling potential U(<$>), we may solve 
Eq. (3.19) for f(t,y). The particular solution depends on a mass eigenvalue m and an 



energy E. We may then use this solution f(t,y) to determine T(t) and X{t) through 
Eqs. ( |3.11 ) and ( |3,12D , What results is the 4d spacetime (described by T(t) and X(t)) on 



which a 4d quantum field (j) with mass m and energy E propagates. We are free to rescale 
f(t,y) to impose T = 1, so that X{t) can then be interpreted as the effective FRW scale 
factor. The crucial result to note is that the scale factor depends on the type of field, its 
coupling potential, and its J^d mass and momentum. 

3.2 Localised fermions 

We now turn to fermions, and perform an analogous calculation to determine the effec- 
tive scale factor describing the 4-dimensional spacetime on which a localised fermion field 
propagates. 

For a 5d fermion fy(t,x l ,y), the action is 

55,* = Jd^xjdy^g- [Vr A E A M (d M + uj m )V - U***] , (3.21) 

where T A are the 5d flat-space gamma matrices, E A M are the vielbeins and com is the spin 
connection 2 . The gamma-matrices obey {r^r^} = 2rj AB , with r\ AB = diag(— 1, 1, 1, 1, 1). 
The coefficient l% of the mass term will in general be a function of other fields, to allow, 
for example, coupling of the fermion to the domain wall. 



2 We are using A, B to denote 5d flat-space indices and M, N to denote 5d curved-space indices. 
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As for a scalar field, we perform a separation of scales in time and separation of 
variables in space 3 

^(t,x i ,y)=u(t,y)^(t,x i ) (3.22) 
and expand the kinetic terms, ignoring it. The action becomes 

S 5 ^ = Jd 4 xj dy na 3 b^ « 2 ^-n^7V + a^^ej djij^j + ... , (3.23) 

where 7 , 7" are the 4d flat-space gamma matrices with {7^,7^} = 2r/ a ^ and ej are 
vielbeins for the 3d space, with 7^ the metric 4 . As before, we require that the correct 
powers of n(t,y), a(t,y) and b(t,y) match with the relevant terms in the prototype 4d 
fermion action 

^7 t0) = / d 4 xT4t)X 3 4t)Vl ? (-TVWV; + X~\th a ejd^) , (3.24) 

where we have used the prototype line element ds\ = —T 2 (t)dt 2 + X^ } {t)^ijdx l dx^ . Match- 
ing kinetic coefficients we obtain 



F^t) = J u 2 a 3 bdy = Xj(t), (3.25) 
GV(i) = I u 2 na 2 bdy = T^(t)X 2 (t) , (3.26) 



which may be inverted to give 



T^t) = F- 2/3 (t)G^(t) =Q u 2 a 3 bdy^j ' (J u 2 na 2 bdy^j , (3.27) 

1/3 



X^(t) = Fl r \t) = ( / u 2 a 3 bdy ) . (3.28) 



V 

These results are similar to the scalar case. As there, we can rescale u(t, y) by a slowly 
varying function of time to enforce XU = 1, so that ip describes a 4d fermion field in a 
spacetime with effective scale factor a^(t) = X{t). Again, the definition of the effective 
scale factor depends on the profile of the particular KK mode tp that one is interested in. 
As before, in the infinitely-thin domain-wall limit u 2 — ► S(by) and T^(t) — ► n(t,y = 0), 
X^{t) — > a(t,y = 0), which recovers the known delta-function brane result. For the 
separable cosmological metric (3.14), is a constant and after normalising u such that 
T^p = 1 we have X^(t) = d(t): the standard result. 

To identify the equation satisfied by u(t,y), we impose the requirement that t/j satis- 
fies the 4d Euler-Lagrange equation with mass and use this to eliminate the spatial 
derivatives of i/j from the 5d Euler-Lagrange equation for ^. This yields 



3 There is a subtlety here: we are assuming that all four components of the Dirac spinor ip have the same 
profile u, which may not be warranted. We expand on this later. 

4 a, /3 — 0, 1, 2, 3 are the 4d fiat-spacetime indices, a — 1, 2, 3 is a 3d flat-space index. 
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where 7 s = i7°7 7 2 7 3 . This has a similar structure to the scalar version ( 3.17 ); in particular 
the tp term quantifies the deviation from ip being a 4d field in the usual, infinitely-thin brane 
definition. 

The appearance of 7 and 7 s in Eq. ( |3.29| ) means localised states on the domain wall 
have an unusual Dirac structure. In the Minkowski-brane case, the 7 cj) term is absent 
and this leads to the localisation of chiral states, which are eigenspinors of 7 s . With the 
presence of 7 , one would naively seek Dirac states which are eigenspinors of both 7 and 
7 5 , which is impossible! It therefore seems that the time-dependent background metric 
leads to unconventionally localised spinor states. To understand this problem more deeply, 
consider seeking solutions to Eq. ( 3.29|) when tp = ipi (and u = ul) is left-chiral, i.e. 



7 ipL = — ipz and m 1 p L = 0. Using a plane wave solution for tpL and expanding its Weyl 
components in order to evaluate 7 ( Vl, we obtain two independent equations for ul 

, ( n' 3a' \ 

L + ( 2^ + 2^ j UL + L = ' (3 ' 30) 

where E^ L is the energy of the chiral plane- wave spinor ipi- 

With a non-trivial background, the only solution to Eq. ( 3.31 ) is ul = 0. Thus, there 
are no localised left-chiral spinors. It may be possible to rectify this problem and find 
localised states which have a certain definite spinor structure by relaxing the separation 
ansatz (recall that each component of the Dirac spinor was assumed to have the same 
profile u, which may be an overly strict assumption), but we will not pursue this line of 
thought further here. 

3.3 The effective Newton's constant 

We shall now briefly discuss how to determine the effective Planck mass, and hence New- 
ton's constant, describing the strength of gravity on the brane 5 . Usually, one expands 
the 5d Ricci scalar in the Einstein-Hilbert action in terms of its 4d counterpart, and the 
numerical pre-factor is identified as the Planck mass. For example, in the RS2 model, one 
uses the metric 

dsj = e~ 2M g$ {x»)dx^dx u + dy 2 , (3.32) 
for which the Einstein-Hilbert action can be written as 

Seh = [ d 4 x [ dy^g M^R (3.33) 



d A x I dy e-^Myf-gW Ml e 2 ^R^ , (3.34) 



D 

where is the 4d Ricci scalar associated with gj$. One then identifies the effective 4d 
Planck mass as 

e -Mv\ dy = ^5 ) (3_ 35) 
-00 A 4 



5 We are concerned here with the definition of the Planck mass for use in cosmological situations, e.g. in 



Eq. (2.17), as opposed to its use in Cavendish- like experiments. 
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which agrees with the result obtained from the effective Friedmann equation in the funda- 
mental brane scenario, Eq. ( p. 13 ). 



Following this approach for domain-wall cosmology, we begin by considering how the 
4d Ricci scalar is embedded in the 5d one. However, there is a problem with this approach, 
because the metric factors for t and x % behave differently at the 5d level, and so the 5d Ricci 
scalar does not separate into a 4d piece plus other terms. To make progress, one might 
consider restricting attention to the 3d Ricci scalar (constructed from the 3d spatial metric 
Jij), which does separate, and identifying its pre- factor in the Einstein-Hilbert action as the 
Planck mass. In other words, we consider just the three spatial components of the metric 
perturbations to determine the gravitational coupling, instead of using the temporal and 
spatial components together. Explicitly, we first write the 5d metric in the general form 

dsl = -n 2 (t,y)dt 2 + a 2 (t,y)£ ij (t,x i )dx i dx j + b 2 (t,y)dy 2 . (3.36) 

The 5d Einstein-Hilbert action can then be expanded as 

Seh = j ^x j dyV^d MlR (3.37) 

D / d A x J dy ncPbyfl M| a^R^ , (3.38) 

where R^ is the Ricci scalar constructed from The goal is to match this to the 4d 
prototype Einstein-Hilbert action associated with the general prototype metric 



dsj = -Ti 4 (t)dt 2 + X l M (t)iij{t, x l )dx l dx J , (3.39) 



which yields 



s(Pf°> = / d'xJ-gW M 2 R^ (3.40) 



D J ' d^xT M (t)X^(t)y/i M 2 X^(t)R® . (3.41) 



By comparing Eqs. ( |3.38|) and ( |3.41 ), we can infer that the 5d theory produces 3d 



(three spatial) metric perturbations with effective Planck mass 

This result requires us to specify the 4d spacetime (by specifying Tm(£) and XmW) before 
we can know the Planck mass. As shown in the previous sub-sections, the 4d spacetime is 
dependent on the particle species, and so we obtain a species dependent Planck mass. This 
may not be so surprising given that each species follows a different line element, but it is 
also possible that the assumption of matching only the 3d Ricci scalar is unwarranted. 

Perhaps a more sophisticated calculation would try to elucidate the effective 4d Ein- 
stein equation, or at least the leading order contribution. Ultimately, we would like to 
identify l/2Mp as the constant of proportionality between the dominant (first order) con- 
tributions to the 4d Einstein tensor t^Gp and the stress-energy tensor ^T^) for a given 
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field in the thin, large-tension brane limit. One possible way to perform this calculation 
would be to analyse the equations of motion for metric perturbations. In the case of 
fundamental-brane cosmology, much of the ground-work for such an analysis has been per- 
formed; see for example [49, 50, 51]. For a domain-wall brane, extra complications arise, 



again due to the averaging of the metric over the extra dimension. Further, it seems that 
in order to identify the 4d metric perturbations, one is forced to perform a separation of 
scales, as in the scalar and fermion case. Such a calculation is beyond the scope of this 



paper, and for our purposes here we adopt Eq. ( 3.42 ) as an approximate definition of the 
effective Planck mass. 

4. Effective scale factor for a thin domain wall 

Having developed the general framework for a domain wall with localised matter fields and 
the associated four-dimensional metric, we would like to better understand the behaviour 
of the effective scale factors and a^. These will, of course, depend on the details 
of the domain wall construction. Furthermore, we need to solve explicitly for the metric 
components n(t, y), a(t, y) and b(t, y) in the presence of this domain wall. We are unable to 
find analytic solutions for a coupled domain-wall gravity system, and numerical solutions 
are beyond the scope of this initial work. To make progress therefore, we will assume 
that the domain wall is extremely thin and that therefore the solutions for the metric 
components are well approximated by the set of equations ( 2.14j ) . 



The profiles of the domain- wall fields will play a role in determining the profiles f(t, y) 
and u(t,y) of the localised scalars and fermions respectively. These localised fields will 
then contribute to the total stress-energy and this back-reaction will modify the metric 
components. However, here we shall ignore such back-reaction effects and consider the 
thin domain wall as a small perturbation to the fundamental-brane scenario presented in 
Sec. ^. In order for this perturbative approach to work, it is necessary that the brane 
localised sources be relatively small, meaning that e<l. 

To compute within this approximation scheme, we first normalise f(t,y) such that 
T = 1, by defining 

f(t,y)=r(t)f(t,y) (4.1) 

so that 

T(t) = r(t) F~ l l\t) &l\t) , X(t) = r(t) F 1/4 (i) G 1 / 4 ^) , (4.2) 



where F(t) and G(t) are defined as in Eq. (|3.7|) but with f(t,y) replaced by f(t,y). En- 
forcing T = 1 gives r(t) = F 1 / 4 (t) G~ 3 / 4 (t) so that X(t) = F x / 2 (t) G -1 /»(t). We may then 
compute F(t) and G(t) by substituting in for the bulk metric solutions (|2.14j ) yielding, for 
example, 

G(t) = Jfandy = a Jf e" 2 ^ 1 - (e + e)e~ M sinh(//|y|) + e e sinh 2 (^|y|) dy . (4.3) 

Requiring that the localisation profile f(t, y) be sharply peaked at the centre of the domain 
wall (y = 0) and fall off rapidly in the bulk translates to f 2 (t, y) sinh 2 (^|y|) — > as y — » ±oo. 
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This condition is consistent with the sufficiently-thin domain- wall brane we are dealing with 
here. Thus, we may ignore the second order term 0{ee) and write 



G(t) = a [h(t)-(e + e)I 2 (t)} 



where 



(4.4) 

(4.5) 
(4.6) 



These integrals, I\(t) and ^2(^)5 are dependent on the exact form of the extra-dimensional 
profile f(t,y). However, if the profile is sufficiently peaked, as we are assuming, we have 
hit) <C Ii(t), because sinh(/i|y|) ~ close to the centre of the domain wall. Under these 
assumptions, we may compute 



h(t) = J P(t,y)e- 2 ^dy , 

h{t)= [ f\t,y)e-^smh^\y\)dy . 



F(t) = al[h(t)-(3e-e)I 2 (t)] , 



so that the effective scale factor for the scalar field becomes 



i(t) = X(t) = oo(t) 
= ao(t) 



1 - 
1 + 



{€ - e) m\ 

e h{t) 



H h(t) 



(4.7) 

(4.8) 
(4.9) 



This is one of the main results of our paper - an explicit, quantitative computation of the 
corrections to the effective 4-dimensional scale factor arising from considering a domain- 
wall brane, rather than a fundamental one. The corrections are proportional to the ratio 
between the rate of change of energy density on the brane and the brane tension, and 
inversely proportional to the effective Hubble parameter. The corrections also depend in a 
non-trivial way on the specific localisation profile of the associated field, so that different 
fields are corrected differently. 

The expression (4.9) satisfies — * ao for the independent limits of a Minkowski brane 
with no sources (e — > 0), and an infinitely-thin brane (I2 — ► 0). For a concrete example of 
this latter limit, consider the profile 6 



fHt,y) = ^^±il^osHwT 2w 



(4.10) 



which obeys f 2 — ► 6(y) as w — > 00 (the thin domain-wall limit) and is a typical example of 
smooth localisation factors (see, for example, El). It is then straightforward to compute 



h(t) 2^T(w + 1) - 4I> + i) 



w 



00 



'TTW 



(4.11) 



6 The time dependence of this sample / would arise from the time dependence of the parameter w, 
corresponding to the brane thickness changing over time. 



- 16 - 



which vanishes in the infinitely-thin wall limit. A better approximation for this quantity 
can be found by solving the differential equation ( p. 19 ) for f(t,y), using the background 
metric components n(t,y) and a(t,y). 

For a fermion field the result for the effective scale factor is almost identical to the 
scalar case 

e J 2 (t)' 



cty(t) = a (t) 
where the relevant integrals are 



1 + 



H Ji(t) 



Ji(t)= / tf{t,y)e-^\y\dy 



■h(t) 



if 



'(t,y)e 2 ^lsinh fj,\y\dy , 



(4.12) 

(4.13) 
(4.14) 



and u(t,y) is defined in a similar way to f(t,y). 

The results from this section, namely Eqs. (|4.9| ) and ( 4.12 ), are concrete expressions 
for modifications to cosmology in a domain- wall brane construction, and are the starting 
point for an analysis of the constraints on such theories from observations. We expect that 
a species-dependent scale factor should have an impact on a vast array of cosmological 
observables, including the predictions of BBN, the era of recombination and the spectra of 
the microwave background and large scale structure. Acceptable cosmological behaviour 
should imply constraints on the brane tension a, which appears in e and /i, and the width 
of the domain wall, which enters implicitly through the localisation profiles / and u. 



5. Conclusions 



A self-contained and self-consistent construction of a five-dimensional theory requires that 
the four-dimensional brane on which the standard model fields are confined be formed as 
a domain wall. Such an object has finite thickness, and couples in a variety of ways to the 
particle physics and gravitational fields of the theory. 

Dimensional reduction from five down to four dimensions involves integrating over the 
extra dimension in the 5-dimensional action. Each term in the resulting 4-dimensional 
effective action then contains numerical factors arising from these integrals over the extra- 
dimensional profiles of the fields. These factors must be absorbed into the fields in order to 
normalise the kinetic terms. In the case of a Minkowski metric on the brane, this absorption 
can be carried out via a single definition, to obtain an effective field theory applicable to 
all fields. 

In this paper we have considered the problem of cosmology on a 4-dimensional domain- 
wall brane. In contrast to the procedure of obtaining a Minkowski metric, this is a compli- 
cated process and, in particular, during the normalisation of the kinetic terms, one must 
take into account the fact that different fields may feel different spacetimes. Thus, for 
domain-wall branes it is not sensible to ask the question "what is the effective scale fac- 
tor on the brane?" . Since the brane has non-trivial dependence on the extra-dimensional 
coordinate y, and since the metric components n(t, y) and a(t, y) are not proportional to 
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each other in the bulk, each 4-dimensional slice at constant y corresponds to a different 
spacetime. The effective 4-dimensional spacetime for a localised field with a smooth profile 
in y will thus be an average over all the different slices. This produces a kind of "dimen- 
sional parallax" effect, since different fields have different averages, and thus a different 
"perspective" of the cosmological evolution of each slice. 

Therefore, rather than seeking the effective scale factor on the brane (which was possi- 
ble in the fundamental brane case), we must instead ask: "what is the effective 4d spacetime 
in which a given localised field propagates?" . For each low-energy 4-dimensional field (each 
species and each mode of the KK tower), we may answer this question by determining the 
effective 4-dimensional line element ds\. If this line element takes the form of an FRW line 
element, then we may define an effective scale factor for the associated field. This is as 
close as we are able to come to answering our original question. 

In the case of a localised scalar field, the effective scale factor is given in general by 
Eqs. ( [3.12 ) and ( |3.13 ), and for fermions one obtains (3.28). For the case of a domain-wall 
brane which is thin enough such that one can approximate the metric components n(t,y), 
a(t, y) and b(t, y) with the solutions from the fundamental-brane scenario, the effective 
scale factors for scalars and fermions are ( |4.9[ ) and ( 4.12; ) respectively. All of these results 
reduce, in the infinitely-thin domain-wall limit, to the results obtained for fundamental 
brane cosmology, where the effective scale factor is the bulk metric component a(t, y) 
evaluated at y = 0. 

Beyond this basic difference between the fundamental-brane and the domain-wall case, 
there are a number of other interesting consequences of our construction. In the cosmo- 
logical scenario, the non-trivial averaging of the metric over y means that, just as different 
modes of a KK tower have different extra-dimensional profiles, it is also true that different 
energies of the same mass have different profiles. Thus, particles with different energies 
feel a different scale factor! An outstanding problem is that of defining a unique Planck 
mass (if possible) at the 4d level. We have provided some insight into the solution to this 
problem, in the form of the initial approximation ( 3.42 ), but a full treatment is beyond the 
scope of this paper. 

The novel features we have presented - species-dependent scale factors and non-unique 
Planck masses - will lead to potentially observable cosmological phenomena. We are cur- 
rently performing a detailed phenomenological analysis to determine how these effects 
constrain parameters such as the width and tension of the domain wall, and the extent to 
which they may allow new approaches to cosmological problems. 
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A. Einstein tensor components 



In this appendix we give explicit expressions for the non-zero components of the Einstein 
tensor 

(A.l) 



Gmn = Rmn — ^9mnR , 



for the case of our metric ansatz, equation ( |2.2[ ). The indices and 5 correspond to t and 
y respectively, while both i and j run over the three spatial dimensions. 
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